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ABSTRACT 

We identify a sample of 22 host galaxies of Type la Supernovae (SNe la) at redshifts 
0.95 < z < 1.8 discovered in the Hubble Space Telescope (HST) observations of the Great 
Observatories Origins Deep Survey (GOODS) fields. We measure the photometry of 
the hosts in Spitzer Space Telescope and ground-based imaging of the GOODS fields to 
provide flux densities from the U -band to 24 jum. We fit the broad-band photometry 
of each host with Simple Stellar Population (SSP) models to estimate the age of the 
stellar population giving rise to the SNIa explosions. We break the well-known age- 
extinction degeneracy in such analyses using the Spitzer 24/j,m data to place upper 
limits on the thermally reprocessed, far-infrared emission from dust. The ages of these 
stellar populations give us an estimate of the delay times between the first epoch of 
star-formation in the galaxies and the explosion of the SNe la. We find a bi- modal 
distribution of delay times ranging from 0.06 - 4.75 Gyrs. We also constrain the first- 
epoch of low mass star formation using these results, showing that stars of mass < 8M Q 
were formed within 3 Gyrs after the Big Bang and possibly by z ~ 6. This argues 
against a truncated stellar initial mass function in high redshift galaxies. 

Subject headings: Galaxies: Evolution, Galaxies: General, Stars: Supernovae 

1. Introduction 

It has been shown that SNe la are standardizable candles in t hat there is a very tight relation 



between the peak luminosity and the width of the light-curve (jPhillipsI Il993l ). For this reason 
SNe la have been used to measure the cosmic distance scale and to probe the expansion properties 
of the universe, leading to the discovery that its exp ansion rate is acce l erating due to the exist ence 



of a repulsive force, given the name Dark Energy (jRiess et al.l Il998l ; iPerlmutter et al.l ll999). It 



is important to fully understand how SNe la are produced in order to minimize the systematic 
uncertainties associated with these measurements and enable the properties of the Dark Energy to 
be investigated as a function of cosmic time. 
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It is widely believed that SNela are the ex plosion of a White D warf (WD) star that has 



reached a critical mass (the Chandrasekhar mass. IChandrasekharlll93ll ). The physical mechanism 
leading to the WD reaching this mass is still highly debated. Broadly, there are two leading ideas 
for SNIa progenitors. The first is that the WD star gains mass via accretion of material from 
a normal companion star which has filled its Roche Lobe, referred to as the Single Degenerate 
(SD) scenario. The second, the Double Degenerate (DD) scenario involves the me rger of two 



WD stars after formation and ejection of a common envelope in a binary system (see iLivio 
Podsiadlowski et al. 2008, for a review). 



2001 



In order to try to differentiate between these possible scenarios it is important to constrain 
the SNIa delay time, i.e. the time between the formation of the stellar system and the supernova 
explosion. Evidence for a range of delay times has existed for some time. For example, observations 
show that SNela are preferentially found in late type rather than early typ e galaxies, suggesting 
they are associated w ith young stellar pop ulations with ages of ~ 50 Myrs (jvan den Berg hl ll990l : 
Mannucci et al.1 120051 ) . IWang et al.1 (|1997l ) showed that SNe la are more likely to be found in 



the disk of a galaxy rather than the bulge, indicating an association with recent star formation 
and short delay times. Studies which calculate the delay time by convolving an assumed star 
formation histor y (SFH) with a delay time distribution (DTD) have suggested longer delay times 
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of a f ew Gyrs (IGal-Yam &: M aoz 



Dahlen et al 



2004: 



Strolger et al. 



20081 : IStroker et al.1 12010 ) . However 



2004 



Dahlen et al 



Forster et al 



2004 



Barris &: Tonrv 



(|2006l ) show that these 



analys es depend strongly on the assumed SFH, giving large systematic errors, prompting lOda et al 



(|2008l ) to fit both the SFH and the DTD. However, those authors were only able to place weak 
const raints on the DTD. Fr om a spectroscopic study of the star- formation histories of local SN la 
hosts iGallagher et al.l (120051 ) put a lower limit on the delay time of 2 Gyrs. 



Mannucci et al.l (|2006l ) have shown that a combination of observations at high and low redshift 



cannot be matched by a DTD with a single delay time, but that they are best matched by a 
bi-m o dal DTD. This hypothesis is developed fur ther by various groups (jScannapieco k, Bildsten 
20051 : 1 Sullivan et al.1 120061 : iNeill et al.1 120061 . 120071 ) who model the SN la rate as a two-component 
distribution with a delayed component dependent on the host galaxy stellar mass and a prompt 
component dependent on the host galaxy star f ormation rate. Several authors have u s ed this model 
to reproduce the observed S N la rates at z ~ 1 (jAubourg et alJl200a ; INeill et alJl2007l ; iDahlen et al 



2004 : Botticella et al. 20081), although again t hese results depend on the assumed SFH. In order to 
avoid such an assumption. iTotani et al.l (|2008l ) used an SED fitting technique to determine the ages 
of the host galaxy stellar populatio ns to d e rive a DTD which is a power-law in the range 0.1 - 10 
Gyrs. From a spectroscopic study. iHowelll (|200ll ) showed that sub-luminous SNela tend to come 
from old populations whereas over-luminous SNe la are from young populations, su ggesting different 
progenitor scenarios for the two populations. Furthermore, iPritchet et al.l (|2008l ) showed that the 
single degenerate scenario alone is not sufficient to explain the observed DTD. By stud ying the local 



enviro nment of SNe la (rather than the properties of the whole of the host galaxy) iRaskin et al 
(J2009J) found that the average delay time of nearby prompt SNe was ~ 0.3 — 0.5 Gyrs, although 
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their model allows ~ 30% of prompt SNela to have delay times shorter than 0.1 Gyrs. 



From a study of the X-ray properties of nearby ellipticals iGilfanov h Bogdanl (|2010i ) suggest 
that the X-ray flux from these galaxies is consistent with only ~5% of SN la arising from accreting 
white dwarfs in elliptical galaxies - such systems are expected to produce X-ray emission for a 
significant time, whereas in the white dwarf merger (i.e. double degenerate) scenario, the X-ray 
emission is only present shortly before the explosion. They do not constrain the progenitors of 
SN la in late-type galaxies where the contribution from the single degenerate scenario could be 
significant. Using a maximum-likelihood inversion procedure. iMaoz et al.l (|2010i ) recover the DTD 
for a sample of local SNe la utilising the SFH of each individual galaxy, finding evidence for both 
prompt (with delay times < 0.42 Gyrs) and delayed (with delay times > 2.4 Gyrs) SNe, where 
the DTD has a peak at sh ort delay times but a broad distribution to longer delay times. Similar 
conclusions are reached by lBrandt et al.1 (120101 ) . 



Greggid (|2010l ) build on the parameterisation of iGreggid (|2005l ) of the evolution of binary 
systems to calculate some features of the DTD of single-degenerate, double degenerate and mixed 
model progenitor scenarios. All models are consistent with both prompt and delayed SNe, although 
the distributions are continuous (suggesting that the distinction is arbitrary). They do explore 
scenarios where there is a mix of SD and DD and where the SD contributes more prompt SNe than 
the DD chann el, however, they see n o theo retical basis for this and c omparisons to the observed 
SN la rates of iGreggio &; Cappellard (|2009l ) and ISullivan et al.l ([200a) do not favour the mixed 



scenarios. 



The Great Observatories Origins Deep Survey (GOODS; iDickinson et al.l [2003) is a multi- 
wavelength survey covering ~ 330 arcmin 2 over two fields. In conjunction with the GOODS survey 
a supernova search wa s conducted, surveying both fields at several different epochs with HST ACS 



Giavalisco et al 



Riess et al 



2004 



2004) resulting in a catalogue of 22 SNe la at z > 0.95 with spectroscopic redshifts, 



2007) . We add a further 3 SNe 1997fg, 1997ff and 2002dd from past SN searches 



in these fields (jGilliland et al.lll999l : iRiess et al 



2001 



Blakeslee et al.1 12003 ). 



In the present work, we compile optical/near-infrared photometry from ACS and NICMOS 
on HST, Spitzer 3.6 — 24 //m data as well as supplementary ground-based data where available for 
the host galaxies of these SNe la. We then fit the multi- wavelength photometry of the hosts to the 
Single Stellar Population (SSP) models of Chariot and Bruzual (priv. com.) in order to find the 
ages of the stellar populations in the SN host galaxies, thus allowing us to study the delay times of 
SNela. 

Furthermore, by calculating the ages of the stellar populations we can constr ain the first epoch 
of low mass star formation (the progenitors of SNe la have masses < 8M0, see e.g. iBlanc Greggio 
2008, and references therein), i.e. the time after the Big Bang that stars must have formed in order 
to yield the stellar populations we find. This allo ws us to constrain mod els which suggest that only 
stars with > IOMq might have formed at z ~ 6 (ITumlinson et al.ll2004l ). 



The remainder of this paper is arranged as follows. In section 2 we describe the identification 
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of host galaxies, in section 3 we outline the SED fitting method. We present our results in section 4 
and provide a discussion of these results in section 5. Finally our conclusions are given in section 6. 
In the rest of this analysis we use the concordance cosmology of $7 a/ = 0.3, = 0.7 and Hq = 70 
km s _1 Mpc -1 . We present photometry in the AB magnitude system unless otherwise stated. 



2. Host Galaxy Identification 



We us e the supernova s ampl e of lRiess et al.l (|2007l ) with additional SNela from iGilliland et al 



(|199a ) and lBlakeslee et all (|2003 ). giving us 25 SNela at z > 0.95. We cross-matched these SNela 
with t he GOODS A CS v2.0 catalogue in both the Northern and Southern Fields using the software 
topcat (jTavlorjbooj l and a cross-matching radius of 1". This resulted in a catalogue of 22 SNe la host 
galaxies. We then match the ACS host galaxy positions to the ground-based and IRAC positions 
with a radius of 0'/5, after correcting the IRAC and ground-based catalogues for the well-known 
/ .'38 offset in the GOODS North catalogues. The SN positions are also corrected for this offset 
where necessary. The large radius used for our initial cross-match ensures that all possible SNela 
with detected ACS host galaxies are identified, however, it could lead to false identifications with 
other galaxies. Therefore we visually inspected the ACS images for each SNIa, but all appeared 
to be good identifications. The three SNe with no matches are 2002dd, HST04Sas and HST04Gre. 
In the cases of 2002dd and HST04Sas, the host is not found in the catalogue due to confusion 
with nearby bright galaxies. HST04Gre appears to be a hostless SNIa, the F-band postage stamp 
is shown in figure Q] along with the host for SN 1997ff in the same band, where the host is well 
identified. The co-ordinates of the SNe la and their host galaxies are given in Table [T] (SNe la with 
ambiguous host identifications are shown with only the SN coordinates given). 



2 r ^ 





L ■■ ■ ■ ■ I 



Fig. 1. — HST V band images centred around two SNela, the positions of which are indicated by 
the red cross. The left hand panel shows the host galaxy for 1997ff is well identified whereas the 
right hand panel shows that the host for HST04Gre has not been identified. It is therefore removed 
from the sample of host galaxies which are analysed in this work. 

In order to break the well-known degeneracy that exists between age and extinction in fitting 
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Table 1. J2000 coordinates of the SNela and the offset to the host galaxies along with the 
spectroscopic redshift of the SNe la (either spectroscopy directly of the SN la or from spectroscopy 
of the host). SNela where the host galaxy offset is missing have a confused or undetected host 

and are therefore dropped from the sample. 
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^Ty pographical error in the co-ordinates given in iRiess et al 



(|2007j ) discovered by comparison to their postage stamps. Orig- 
inal co-ordinates are HST04Eag 189.33646, 62.22820; HST04Mcg 
53.04175, -27.83055, corrected co-ordinates are those presented here. 
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model spectral energy distributions to optical/NIR photometry, we use the MIPS 24^m photometry 
to place constraints on the fraction of energy that is absorbed by dust and re-emitted at longer 
wavelengths. This results in upper limits on the amount of extinction applied to the model SED. 

For each host we then have HST ACS BViz, Spitzer IRAC chl - ^ and MIPS 24/mi 1 data as 
well a s ground basec@ UJHK. We also have some HST NICMOS JH coverage from Buitrago et al 



(|2008l ). We reject photometry with statistical errors larger than 0.3 mag which are indicative of 
marginal detections. We then add a systematic err or in quadrature to the photometric errors given 



in the catalogue of 5% for HST ACS and NICMOS (|Pavlovskv et alJl2004l : lThatte et al 
for IRAC chl — 3 and ground-based photometry and 15% for IRAC ch4 (jReach et al 



2009), 10% 



20051 ). We 



remove the southern U band data from the analysis due to inhomogeneous variations across the 
field. 



The process of mosaicing the ACS data includes the flux from the SNe in the lRiess et al.l ([20071 ) 
sample. We have estimated and subtracted out the SN contribution to the host galaxy photometry 
in th e catalogue using the following process. We use the SN photometry given in lRiess et al.l (|2004l . 
20071 ). In each band we add up the total SN flux and divide by the number of epochs of observations 
going into the ACS stack and subtract this from the host flux. This effect is present in the i and z 
band s, with a small con tribution in the V-band (partly due to the SN colour and partly because the 
later iRiess et al.1 (|2007l ) SN search did not include F-band re-imaging of the GOODS fields). The 
5-band data was taken at an earlier time to the Viz data and is thus not affected. In several cases 



Riess et al.l (120041 ) give i and z band SN photometry but not the F-band photometry. In these 
cases, we estimate the contribution to the V— band by extrapolating a simple power-law spectrum 
fit to the nearest two bands. In a very few cases, only z— band photometry is given; in t h ese ca ses 
we fit the power-law to the observed colours of Type la SNe as presented in |jha et al.l (|2007l ) to 
measure the contribution of the SN in the i— band and the V— band. In all cases the correction to 
the host photometry in the F-band is small (there are only 3 cases where the correction is larger 
than 0.01/iJy with the largest being 0.076^Jy). In the i and z-bands the correction is small in most 
cases, the correction is larger than ~ 10% of the host flux in only 4 and 5 cases respectively. In the 
cases of HST05Gab and 2002fw the SN contamination is very large, giving unreliable photometry 
in the i and z bands. Removing these bands leaves us with good photometry only in the BV, IRAC 
chl and ch2 in the case of HST05Gab and only BV and K in the case of 2002fw. Since we are 
unlikely to be able to constrain the host SED with only 3 or 4 bands we remove these SNe from 
the analysis. In the remaining cases a visual inspection further suggests that the contamination is 
small. Therefore, although the details of the stack process could affect the precise correction that 



1 with the exception of SNe HST04Pat, which was off the edge of the complete Spitzer survey, and 2002ki which 
had an uncertain IRAC counterpart due to confusion 

2 Ground-based data is from ISAAC/Very Large Telescope, Keck Telescope (spectroscopy), MOSAIC/Kitt Peak 
National Observatory, WIRCAM on the Canada-France Hawaii Telescope and the MOSAIC/Cerro Tololo Inter- 
American Observatory. The full list of data sets can be found at http: / /www.stsci.edu/science/goods / 1 
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is required, this effect is likely to be small. 



The final photometry used, including the systematic errors added in quadrature is given in 
tables [2] (HST ACS, ground-based) and [3] (Spitzer IRAC and MIPS). The photometry used i n 
the present work for the host g alaxy of SN 1997ff differs from that given in iRiess et al.l (|200ll ). 
In particular, IRiess et al.1 (|200ll ) give i?-band photometry of 26.67 ± 0.16 mag, whereas we do not 
find an accurate detection for the host-galaxy in this band (the catalogue value is 28.23 ± 1.26 
mag so it is remove d from our anal y sis du e to the large error). Further investigation revealed that 
the host is given in I Williams et al.l (| 19961 ) as 5-band magnitude of 29.18 with a signal-to- noise of 
3.5. A similar pattern is seen in the iz bands suggesting that perhaps there is a variable Active 
Galactic Nucleus (AGN) in this host, although further investigation beyond the scope of this work 
is required to confirm this. 



Table 2. HST ACS and ground-based photometry of the host galaxies used in the fits, the errors given include additional 

systematic errors (see text). All magnitudes are in AB magnitudes. 



SN Name U [mag] B [mag] V [mag] i [mag] z [mag] J [mag] H [mag] K [mag] 
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Table 3. Spitzer IRAC and MIPS photometry of the host galaxies used in the fits, the errors given include additional 

systematic errors (see text). 
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13.4± 1.3 


9.0± 0.9 


6.4± 0.7 


3.4± 0.6 


< 25 


HST04Eag 


10.3± 1.0 


7.6± 0.8 


5.8± 0.7 


5.7± 0.9 


50.9± 5.1 


HST05Gab 


0.6± 0.1 


0.3± 0.1 


< 1.5 


< 1.5 


< 25 


HST05Str 


5.4± 0.5 


3.8± 0.4 


2.8± 0.5 


2.3± 0.6 


< 25 


1997fg 


7.8± 0.8 


5.6± 0.6 


4.5± 0.5 


3.7± 0.6 


< 25 


2003aj 


1.6± 0.2 


1.3± 0.1 




< 1.5 


< 25 


2002fx 


0.4± 0.1 


0.4± 0.1 


< 1.5 




< 25 


2003ak 


5.0± 0.5 


5.1± 0.5 


3.3± 0.5 


2.2± 0.6 


< 25 


2002hp 


35.4± 3.5 


31.9± 3.2 


22.2± 2.2 


15.4± 2.3 


< 25 


HST04Mcg 


20.8± 2.1 


20. 1± 2.0 


14.7± 1.5 


13. 1± 2.0 


59. 5± 6.0 


HST04Omb 


4.3± 0.4 


3.0± 0.3 


1.8± 0.4 




< 25 
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3. Model Fitting 



We use the population synthesis mod els of Chariot and B r uzual (CB07, priv. com.) which are 
generated from an updated version of the iBruzual &; Charlotl (|2003l ) GALAXEV code to produce 
models which include a new prescription for Thermally Pulsating AGB stars. The code calculates 
the spectral evolut i on of a stellar population based on a library of observed stellar spectra (see 
Bruzual fc Charlotl ([2003) and references therein). We generated a suite of models with varying 
metallicity, Initial Mass Function (IMF; i.e. the distribution of stellar masses for the starburst), 
stellar population age and star formation history (SFH, either a constant SFR or an exponentially 
declining SFR wit h e- folding tim e r) . For each model we apply a dust extinction correction of 
varying Ay with a lCalzetti et al.l ( 2000 ) dust extinction law, giving a total of ~ 5 x 10 5 models. 
Table U] shows the parameter values used to generate the models. Finally, we ensure that the age 
of the stellar population never exceeds the age of the Universe at the redshift of the SN la. 



Table 4. Input parameters used to generate the model SEDs 



Parameter 


Allowed Values 


Metallicity 


0.005, 0.020, 0.200, 0.400, 1.000, 2.500 Z Q 


IMF 


Salpeter: ^ = M~ 135 O.1M < M < 100M o 

Chabrier: ^ = exp(-(log 10 (M) - log 10 (0.08)) 2 /0.9522) 0.1M Q < M < 1M Q 
= M- 1 - 3 O.1M < M < 1OOM 


SFH 


ijj(t) = r- 1 exp(-i/r), r = 0.01, 0.03, 0.1, 0.3 1 Gyrs and r = oo i.e. constant SFR 


Stellar population age 


unevenly spaced in the range 0.001 - 6 Gyrs 


Dust, Ay 


0-5 mae. step size of 0.1 mag. starburst extinction law CCalzetti et al. 2000) 
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Each model is then converted to the observed frame using the spectroscopic redshift of the 
SNIa and convolved with the response/filter functions of the instruments in order to generate the 
equivalent observed photometry for each model in each band. We then calculate the x 2 °f each 
model compared to the data, according to the equation 



* 2 = £^? 



a? 



where f° is the observed flux in each band, f m is the model flux in each band, a is the observed 
error (including the additional systematic error) , b is a normalization factor (calculated as the mean 
ratio of the observed flux to the model flux, weighted by the errors) and the summation is over all 
bands. We then find the model with the minimum \ 2 for each host galaxy. 

In order to break the degeneracy between young models with large extinction and old models 
with low extinction we use the 24/im photometry since it is a reliable proxy for the fraction of 



light that is absorbed by dust and reprocessed at longer wavelengths. By using the lCharv &: Elbaz 



(|200ll ) templates, the 24^m photometry can be translated to a total far-infrared (FIR) luminosity. 
This provides an upper limit to the fraction of optical/near- infrared light that is absorbed by dust. 
We then compare this to the FIR luminosity of each model, calculated thus 

1 = J l x (l - exp(- J AvK/1.086))dA 

where l\ is the flux o f the model at each w avelength A and n is the extinction correction given by 



the extinction law of lCalzetti et al.l ((200fJ). Any models which have an FIR luminosity that is in 
excess of the FIR luminosity calculated from the 24/im photometry /limits are rejected. We use 
upper limits to reject any models that exceed the photometric limits. The limits we use are 28.6, 
28.6, 27.9, 27.4 AB mag for the ACS BViz bands (5a limits), 0.5, 0.5, 1.5, 1.5/jJy for the IRAC 
chl-4 and 25/uJy for the MIPS 24/im data (80% completeness limits). 

Each SED in the CB07 models is the combination of the SEDs from all the various stars that 
have formed (and evolved) over the lifetime of the model galaxy. All but one of our star-formation 
histories is exponential, in these cases most of the stars are formed at the redshift of formation. 
The final time-step from the models is the upper limit to the age of a star in the population. While 
most of the stars in the population are old, most of the luminosity comes from younger stars. In 
order to account for this we calculate a mean age weighted by the fractional contribution to the 
y-band luminosity from the stars of each age. 



4. Results 

4.1. Delay Time Distribution 

We have found the best-fit stellar population model for each host using a minimum x 2 tech- 
nique. The set of best fitting parameters for each host along with the x 2 °f the fit is shown in 
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Table [5j Since we have calculated the x 2 f° r each model, we calculate the errors by finding those 
bins with x 2 < Xmin + Xgaussi where Xmin is the minimum x 2 of the fit and x 2 ga uss is tne X 2 value 
for a 68% and 95% confidence level for a normal x 2 distribution with the same number of degrees 
of freedom as used in the fit (given by the number of photometric points less the number of pa- 
rameters). The best-fitting SEDs along with the observed photometry used in the fit are shown in 
figure [21 It is also worth noting the effect of using the 24 //in limit at this stage. Figure [3] is an 
enlarged version of the SED plot for 1997fg but also shown is the best-fitting SED found if the FIR 
luminosity limit is not used. The figure shows that both SEDs are reasonable fits to the data (with 
X 2 of 11.1 when the limit is included and 3.6 when it is not), however by including the limit, the 
best-fit luminosity-weighted age changes from ~ 0.02 Gyrs to ~ 0.09 Gyrs. 



Table 5. Table showing the best-fit parameters for each SN host (labeled by SN). Errors are calculated by finding the 
parameter value which raises the measured x 2 (marginalised over all other parameters) above the 68% confidence threshold as 
calculated for figure [6l In cases where the best-fit is at the extreme end of parameter space or if the marginalised x 2 never 
reaches larger than the threshold no error is given. (Age), Z and M are the luminosity weighted age, metallicity and stellar 

mass respectively 



SNIa 


A 


.ge [Gyrs] 


(Age) [Gyrs] 


Ay [mags] 




Z[Z ] 


r[Gyrs] 


IMF 


xl 


M [logio(M Q )] 


1997ff 





+0.444 
ol 1 -0.062 





r fi1 +0.341 
Ju -0.090 


5 +0 - 4 

U -°-0.3 


2. 


500-1.500 


001+0-09 


salpeter 


3.05 


11.06lg;«? 


2003es 


2 


OQQ + 0.700 

Juu -1.161 


1 


QQQ+0.700 

• 3a3 -0.898 


3+ 47 

u -°-0.3 


1 


Ooo+i-soo 

uuu -0.995 


0.30-0.29 


chabrier 


1.95 


lo.sotS 


2003az 


4 


750-3.316 


3 


oq 7 +0.001 
• 8y ' -2.763 


n 7 +0.6 
u -'-0.1 


2. 


500-2.100 


1.00-0.70 


salpeter 


4.71 


11 i 9 +0.301 
n.i__ 673 


2003dy 


0. 


9Qf;+4.214 
0.172 





18n +l.568 


n ,1+0 2 

0-4Io.4 


0. 


A r\r\-\-2 100 


0.10-0.09 


chabrier 


0.65 


a n i +0 731 


2002ki 





nn i:+4.995 
uud -0.000 





002 +1.961 
" uz -0.000 


1 6+ 34 





005+ 2 ' 495 


0.01 


chabrier 


3.40 


a qo+1.291 
°- ao -0.070 


HST04Pat 





nfi4 +0.017 
UD *-0.014 





054+0-017 
U04 -0.014 


1 «+°- 2 


0. 


005+ 2495 


0.01 


chabrier 


5.22 


1 7q+ 288 

1U. I »_o.071 


HST05Fer 


2 


r m +3.000 
OUU_2 214 


1 


7n c+0.484 
' U8 -1.528 


0.0+ 07 





400+2 100 

/±uu -0.200 


1.00-0.90 


chabrier 


1.83 


Q O o+0.354 
y - zo -0.381 


HST05Koe 





n7o +0.089 
u ' z -0.008 





nfi9 +0.062 
uoz -0.008 


1 S+ 3 - 2 
1 - 8 -0.2 


2 


.500-1.500 


OQl+0.02 


salpeter 


2.02 


10-23-Q 332 


HST05Red 





nc:c:+0.586 
uod -0.023 





no, +0.205 
' uo -0.009 




1 


Ooo+i-soo 

uuu -0.995 


0.03-0.02 


chabrier 


2.55 


8.84l»;« 2 


HST05Lan 


1 


7nn +0.200 
' ' uu -0.091 


1 


40f) +0.200 
•^""-0.091 




1 


Ooo+i-soo 

uuu -0.600 


0.30-0.29 


salpeter 


1.66 


10.661°;°" 


HST04Tha 


2 


onn+3.550 

zuu -1.061 


1 


oqq+2.997 
oyy_ .798 


00 +o.5 


1. 


Ooo+i-soo 

uuu -0.800 


O0+0.70 
U.OU_ 29 


chabrier 


2.02 


10.071°;^ 


HST04Eag 


5. 


.500-5.436 


2 


1Q9 +3.298 
.1»Z_2 138 


o.6l° ;„ 





005+ 2 ' 495 


Constant SFR 


salpeter 


2.49 


10-48lg; 47 ? 


HST05Str 





n79 +0.072 
u ' z -0.008 





nfi2 +0.045 
uoz -0.008 


0-615:! 


2. 


.500-i.soo 


0i+o.02 


chabrier 


1.53 


n 9 r+0.353 
y - zd -0.064 


1997fg 





1n9 +2.898 
±uz -0.102 





aqo+2.088 
•" az -0.092 


4+0.1 

U -*-0.3 


2. 


.500-2.495 


QOl+0.99 


chabrier 


1.86 


O 40+0. 860 
y -*°-0.086 


2003aj 





lfil +4.339 
' ±UJ --0.108 





i 24 +1.624 
lz ^-0.082 


o i+ - 9 


1. 


Ooo+i-soo 

uuu -0.995 


0.03-0.02 


salpeter 


0.77 


q 25+0-611 

y - zd -0.335 


2002fx 


2 


7 r n +1.500 
• ,ou -2.746 


1 


O r,r + 0.419 
ZZJ -1.224 


0.0+ 1 ' 2 





005+ 2 ' 495 


Constant SFR 


chabrier 


1.63 


Q4+0.373 


2003ak 





19 o+0.127 
1ZO -0.078 





1 1 c+0.035 
llo -0.078 


OQ+0.6 


1. 


Ooo+i-soo 

uuu -0.995 


OQl+0.09 


chabrier 


0.53 


q 7n +0.435 

a - ' u -0.160 


2002hp 


1 


7nn +2.800 
' ' uu -0.795 


1 


rn 7 +2.903 
oot -0.720 


3+ - 5 


1. 


Ooo+i-soo 

uuu -0.800 


-,0+0.20 
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salpeter 


0.45 
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ii.uu_ 459 


HST04Mcg 
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salpeter 


0.56 
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Fig. 2. — SED plots for the best fit model for the SNe host galaxies, blue points are optical data, 
green points are near-infrared and red points are IRAC chl — 4. 




figure continued 
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Since we are only fitting one population to the photometry of the host galaxies it is possible 
that the SNIa progenitors were born in a more recent starburst which only contributes a fraction 
to the total emission of the galaxy which is dominated by an older population. In this case we 
would overestimate the ages as we would fit to the dominant, older population but the SNIa in fact 
comes from a younger population. In order to test this possibility we perform a two-component fit. 
We use the best-fitting parameters for metallicity, SFH and IMF and take a model which is as old 
as the universe at the host redshift and a model which is 10 Myrs old. We then add the two models 
together allowing the fractional contribution from the young population to vary between - 1 in 
steps of 0.1. We also allow the extinction to vary and to be different in the two populations. In all 
cases the resulting fit is worse than the original best-fit using only one population. This suggests 
that our assumption that the SNe la most likely come from the one population model we are fitting 
is valid. 



We also assess ed the ICharlot Fall (|2000l ) prescription for the dust ext inction. We find tha t 
the ages using the ICharlot k, Falll (|2000j) recipe differ from those using the ICalzetti et al.1 (J2000J) 
recipe by more than a factor o f 2, in only 3 cases. In two of these cases these ages are larger than 
those u sing the ICalzetti et al.1 (j2000l ) recipe, however, in all 3 cases the x 2 of the best-fit is lower 



using a 



Calzetti et al. 



(2000) recipe compared to the ICharlot h Falll (|2000l ) recipe. We conclude 
that our derived ages are robust with respect to the fitting technique. 

Figure [4] shows contour plots of the x 2 distribution in the age - Ay plane for the SNe host 
galaxies, showing 68% and 95% confidence intervals. The figure shows that in some cases there is 
a fairly smooth distribution of the degeneracy between the age of the stellar population and the 
dust extinction, in that good fits can be achieved with younger populations with a higher dust 
extinction, as expected. This degeneracy remains despite our removing any models which give a 
FIR luminosity that is inconsistent with the MIPS 24/mi data, although it is much reduced. Many 
panels of figure [U however, do not show a smooth distribution. While they generally show the same 
degeneracy, it is clear that the age parameter is not sampled sufficiently, in particular at older ages 
to give a smooth distribution. Unfortunately, we are unable to alter the ages of the output of the 
stellar population models. However, the uncertainties in tableland figures (U [10] span the entire 
range of ages that give consistent fits. 

Figure [5] shows the luminosity weighted stellar population ages of the host galaxies. This age 
is plotted for each SNIa host with the reduced x 2 (x 2 = X 2 / v i where v is the number of degrees of 
freedom) of the best-fit model. The error bars are calculated from the grid of x 2 > whereby we find 
the age at which the x 2 rises above the 68% confidence level, collapsing over all other parameters. 
The plot shows that there are two hosts with exceptionally large error bars, namely the host galaxies 
of 2002fx and 2002ki. In the latter case this is most likely due to the lack of IRAC photometry, in 
the former case this is most likely due to large photometric uncertainties. In any case we remove 
these hosts from the analysis (including the calculations based on this plot). The plot shows a large 
range of luminosity weighted ages from 0.03 - 3.90 Gyrs, with both prompt and delayed SNe la. 
As we wish to constrain the maximum delay time, we show the plot again but when the upper-age 
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100 f 




0.1 1.0 10.0 



Fig. 3. — SED plot for the best fit model for the host galaxy of SN la 1997fg, blue points are optical 
data, green points are near infrared and red points are IRAC chl — 4. The black line shows the 
best-fit model SED, the magenta line shows the best-fit model SED when the FIR luminosity limit 
is not included. The best-fit luminosity- weighted age, Ay combination changes from (0.02 Gyr, 2.1 
mag) to (0.09 Gyr, 0.4 mag) when the limit is included. 
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limit (i.e. the best fitting time-step of the CB07 models) rather than the luminosity-weighted mean 
age is used in figure [6) Both plots show that the weighted mean age of the SNe hosts is < 0.1 Gyr 



i mply ing that they harbor young stellar populations. This fact was noted earlier by I Chary et al 



(|2005l ) who found that ~60% of the SN la host galaxies harbor young, dusty starbursts based on 



their detection at 24 /im; this is a factor of 1.5 higher than the field galaxy population. 

The distribution of the best-fit luminosity-weighted stellar population ages are shown in fig- 
ure [7J We show the his togram when only high confidence SNe la are included, these are SNe which 



are classed as 'Gold' in lRiess et al.l (|2007l ) and also have good spectra of the SNe (red-dashed line). 
The black solid line represents all SNe and includes 'Silver' SNe and 'Gold' SNe without spectra. 
'Bronze' SNe are not included in any of this analysis. The figure suggests a bi-modal distribution 
with median ages of 0.11 and 1.9 Gyrs for the two populations (both the whole sample and when 
only considering the high-confidence SNe la). We arbitrarily consider a host age < 0.4 Gyrs to be 
young, as this is the minimum in the distribution and since there is no precise definition in the 
literature. We perform a KS-test to determine whether a bi-modal distribution is a good descrip- 
tion of the data. Comparing to a distribution comprised of two gaussian distributions centered at 
the young and old median ages (allowing the standard deviation of each gaussian and the ratio 
between their amplitudes to vary) gives a KS-statistic of 0.99 for the high confidence SNe la and 
0.97 for the whole sample. The KS statistic for the distribution to be fit by a single Gaussian has 
a probability of only 42%. Thus, a bimodal distribution is clearly favored by the age distribution 
although a larger sample of objects such as those available from the upcoming CANDELS/WFC3 
survey is required for definitive confirmation of a bimodal distribution of ages. 

Due to the large error bars for the ages of the stellar populations we perform a monte carlo 
simulation to test the strength of any bi-modality. For each simulation, we calculate an age for each 
SN from the probability distribution calculated from the x 2 distribution of the error bars shown 
in figure [5j We then repeat the KS-test. We perform this simulation 1000 times. We find that 
a bi-modal distribution is preferred 95% and 97% of the time for the whole sample and for the 
high confidence sample respectively. Finally, we use the simulation to test the existence of prompt 
SNe given the large error bars of figure [H The mean number of young hosts in the simulation 
is 9.5 ± 0.05 and of old hosts is 10.5 ± 0.05 for the full sample. When considering only the high 
confidence SNe la the average number of young and old hosts is 4.6 ± 0.03 and 6.4 ± 0.03. These 
results suggest that we can be fairly confident of the existence of prompt Type la SNe at z > 1. 

Finally, we wish to calculate the delay time distribution (DTD). In order to do this we must 
account for the selection efficiency of the supernova search. Typically, this is achieved in the form 
of a 'control time' calculation (which is the total time that a SN co u ld ha ve been detected). We 



use the control time calculation for this SNe search of IDahlen et al.1 (120081) as an estimate of the 
probability of detecting a supernova as a function of redshift. IDahlen et al.1 (|2008l ) give supernova 
rates and the number of supernovae in four redshift bins between 0.2 and 1.8. We use the three 
highest redshift bins to interpolate the control time at each of the redshifts of our SNe using the 
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equation 



R 



N 
t r AV 



where R is the SNela rate, N is the number of SNe observed in a redshift bin, t c is the control 
time and AV is the volume of the redshift bin (IStrolger et al.1 12010 ). We then divide the best-fit 
luminosity-weighted stellar population ages by this control time. In order to calculate the DTD 
we then calculate a histogram of this distribution and divide the histogram by the equivalent 
distribution of a sample of the field population in the GOODS survey with spectroscopic redshifts 
and z > 0.95. The field sample we use consists of 1507 galaxies across the GOODS-N and GOODS- 
S fields. Figure [8] shows a histogram of these ages and figure [9] shows the resulting DTD. We again 
perform a monte carlo simulation of the DTD. We find that a bi-modal distribution is preferred 80% 
of the time for the whole sample and 74% of the time when considering only the high confidence 
SNela . For the whole sample an exponential distribution is preferred 11% of the time and a single 
gaussian distribution 8.7%. For the high-confidence sample a single gaussian is prefered in 14% of 
the simulati ons and an exponen tial distribution 7% of the time. In figure [9] we also plot the DTD 
obtained by iTotani et al.1 (120081) . For most bin s the two measurements are consistent, however, 
between 0.4 and 2 Gyrs the lTotani et al.1 (|2008l ) resu lt has a large r SNe rate than that found here 
and their results do not appear to be bi-modal. The lTotani et al.l ([20081 ) result was obtained over 
a lower redshift window, extending between 0.4 < z < 1.2 and it is possible the difference reflects 
a change in the dominant SNe la progenitor. 



4.2. Evidence Against Truncation of the High Redshift Stellar IMF 

Since SNela are believed to be the explosions of White Dwarfs (WD) which have reached the 
Chandrasekhar mass, the stars eventually giving rise to these explosions are low mass stars. A star 



of > 8 Mq will explode as a core-collapse supernova instead of forming a WD (jBlanc fc Greggio 



2008). Therefore, with the ages of the underlying stellar populations which produce the SNela we 
can constrain the first epoch of low mass star formation, T s j 

T.sf — T z T s i e n ar 

where T z is the age of the universe at the redshift of the host galaxy and T ste Uar is the best-fit 
luminosity-weighted stellar population age. T s f is the time since the Big Bang when the stellar 
population formed. The results are shown in figure [10] where the error bars are calculated as above. 
The figure suggests that < 8M Q stars formed within 3 Gyrs of the Bi g Bang and possibly by 



z ~ 5. If these results are confirmed they are in contrast to the proposal by lTumlinson et al.1 ([2004) 
who suggest that instead of requiring the first stars to be very massive stars (M > 140Mq) the 
primordial IMF may be truncated at ~ 10 — 20M Q at z > 6. This suggestion was primarily based 
on the Fe-peak and r— process elemental abundance patterns of extremely metal poor stars in the 
Galactic halo. 
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Interestingly, we do not see any trend of peak absolute magnitude with host age; based on the 
pr evious work we would expect younger galaxies to host fainter SNe la . However, as pointed out 
by I Sullivan et al.l (|2010l ) the sample used here is perhaps too small to see the slight trend found by 
those authors. Furthermore, iNeill et al.l (|2009l ). found that for a sample of local SNe la that there 
was no trend with light-curve corrected peak B-band absolute magnitude and host luminosity 
weighted age except when a subset of hosts with low dust extinction was considered. 



5. Discussion 

We have utilised the multi-wavelength spectral energy distribution of the host galaxies of Type 
la SNe at z > 1 to identify a possible bi-modal distribution of luminosity-weighted stellar ages and 
thereby delay times between the burst of star-formation and the time at which the SN explodes. 
We find evidence for both prompt (i.e. short delay times, < 0.4 Gyrs) and delayed (i.e. long delay 
times) SNe la, with some extremely young (< O.lGyrs) luminosity weighted ages. We discuss our 
results in the context of past measurements at low and high redshift. 



5.1. Low Redshift 



In the low redshift universe there is evidence for both prompt and delayed SNe la. lAubourg et al 
( 20081) find sign i ficant evidence for a population of SNela with progenitor lifetimes of <0.18 Gyrs. 



Mannucci et al 



(2005) find that SNela are more common in blue rather than red galaxies and 



Delia Valle et al.l (|2005l ) find that the SN la rate is 4 times higher in radio-loud rather than radio- 
quiet galaxies, sugge sting that SNe la a re associated with younger stellar populations and therefore 
shorter delay times. ISchawinskil (|2009l ) fit SDSS and GALEX photometry of the host galaxies of 
21 local SNela in early type galaxies to a two-component Stellar Population model, with an old 
component (of age varying 1-15 Gyrs) to represent the older, underlying population and a young 
starburst component with varying age and mass fraction. They find no SNe la with delay times 
< 0.1 Gyrs, and a range of minimum delay times of 0.275 - 1.25 Gyrs. This is perhaps to be expected 
as only SNe la in early type galaxies are studied. Furthermore, measurements of the SN la rate at 
different redshifts have sugge sted that the delayed SNela give a more s ignificant contribution to 



the total rate at low redshift (Sul livan et al. 



2006 



Neill et al. 



2006 



20071 ) . 



Gallagher et al.1 (|2005l ) perform a spectroscopic study of the host galaxies of 57 local SNe la to 



deduce the SFR and SFH of the host galaxy. By measuring the Scalo-b parameter they see evidence 
for a bi-modal distribution which further suggests two progenitor classes for the SNe. However they 
put a lower limit on the delay time of 2 Gyrs. It is hard to see how to reconcile these two results, 
although again these are local galaxies when we would expect more delayed SNe than prompt SNe. 



Neill et al.l (|2009l ) perform an SED fitting analysis of UV a nd optical pho t ometr y of the host 



galaxies of a sample of local SNela. They confirm the results of I Sullivan et al.l (|2006l ) who showed 
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that brighter SNe occur in galaxies with higher specific SFR. lGallagher et al.l (|2008) obtained optical 
spectra for the host galaxies of 29 SNe la selected to be local early type galaxies. From comparisons 
to stellar population synthesis models they find a correlation between age or metallicity with peak 
SNIa V— band absolute magnitude, preferring a trend with age (based on the trend of SN Ia rate 
with specific SFR) such that SNe la from older populations are fainter. iHowell et al.l (|2009l ) find a 
very weak correlation between luminosity-weighted age of the host and 56 Ni mass derived from the 
integrated luminosity of the SN also suggesting that older, low mass progenitors produce fainter 
SNe la. 



5.2. High Redshift 



Many authors have attempted to constrain the delay time distribution (DTD) of SNe la by 
comparing t he observed SNIa rate to that predicted by a convolution of the DTD with an as 



Dahlen et al 



sumed SFH, (Gal- Yam & Maoz 2004: Stroker et al 



2004 



Dahlen et al. 



2004 



Barris &: Tonrv 



2006 



20081 ; IStrolger et al.l l2010l ) . These studies have argued for a range of characteristic 



delay times, spanning 1-4 Gyrs. 

This seems at odds with our results which show a large proportion of SNe la with delay times 
< 0.1 Gyrs. Indeed, the IStrolger et al.l (|2004l . l2010l ) result, finding a characteristic delay time of 
3 — 4 Gyrs, is based on an anal ysis of the same set of SNe used in the present work. iDahlen et al 



(120041 . I2008T ) use the models of IStrolger et al.1 (j2004j) with a measurement of the GOODS SNIa 
rate to show t hat th e best fit DTD is a Gaus sian with a mean delay time of 3.4 Gyrs. However, 
Forster et al.1 (120061 ): iBlanc &: Greggid (|2008l ) have shown that the results of such analyses are 
strongly dependent on the assumed SFH, introducing syste matic errors , and as such our results 
are not necessarily at odds with those authors. Results from lOda et al.l (|2008l ) which attempt to 
fit both the SFH and DTD s i multa neously are only able to put weak constraints on the DTD. 
Furthermore, iPoznanski et al.l (120071 ) measure the SNIa rate at a similar redshift range using a 
dataset from the Subaru Deep Field a nd find a more con s tant r ate which could su g gest s horter 
delay times. In a companion paper to iGal-Yam &: Maoa (|2004l ). iMaoz Gal- Yarn! (|2004l ) show 
that iron abundances in clusters require delay times of < 2 Gyrs. 

These seemingly contradictory results, with evidence for both prompt and delayed SNe la a t 
high and low redshift have led to the suggestion of a two-component DTD (IMannucci et al.ll2005l ). 
For example, IMannucci et al.l (|200 6) use several datasets to show that the observations cannot 
be simultaneously matched by a single delay-time. Several authors have developed this further 
suggesting a model with a prompt component dependent on the specific SFR (SFR per unit stellar 
mass) of the host galaxy and a delayed component dependent on the s t ellar mass of the host galax y 
(|Scannapieco fc BildsterJ bood : ISullivan et al.ll2006l : llSTeill et aklbood . I2OO7I : kubourg et al.l Eooil ^. 
These DTDs are then matched to the observed SNIa rate to fit the model parameters, tending to 
find a best-fit DTD dominated by the prompt component, especially at high-redshift. Furthermore, 
in many of these models the contribution of the prompt component is expected to increase with 
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redshift ([Sullivan et al.l 120061 ) . The average SNe la light curve width appears to increase with 
redsh ift, supporting the s e models (prompt SNe la are more luminous and have a broader light 



curve 



Howell et al 



20071 ; ISullivan et al.ll2009l ). However, some studies have shown that SN la rate 



measurements are unable to differentiate between DTD models to any significance (jNeill et al 



20061 . 120071 : Blanc fc Greggicj |200J ; Botticella et al 



ar e still highly depend ent on the choice of SFH. 



2008 



Kuznetsova et al 



Oda et al.l 120081 ) and that these results 



(2008) re-computed the results 



of lDahlen et al.l (j2004l ) using a more sophisticated technique and additional data to find that they 
could not discriminate between a two-component model and a gaussian single delay time DTD. 



Totani et al.l (|2008l ) perform an analysis similar to that presented here using SN la host galaxies 
in the Subaru/XMM-Newton Deep Survey (SXDS) in the redshift range 0.4 — 1.2. They find a 
DTD in the range 0.1 — 8 Gyrs (extending to 10 Gyrs from the SNIa rate in local ellipticals). This 
analysis selects old, passively evolving galaxies and therefore does not probe the shortest delay 
times. However, they do not find a bi-modal distribution as we find here. 

There are a number of implications of this result for cosmological studies. That there is a 
large population of SNe la which have short delay times means that these objects could be used as 
cosmological probes to very high redshift (z > 3). Furthermore, if prompt and delayed SNe la have 
different light-curve shape/luminosity relations the host galaxy will need to be taken into account 
when using SNe la to determine the equation of state parameter of dark energy, w, especially when 
considering p ossible evolutio n of th is parameter. Using a large dataset of SNe la spanning a wide 
redshift range ISullivan et al.1 (|2010l ) showed that SNe la light-curve widths depend on host galaxy 
specific star formation rate (SSFR) and stellar mass, with narrow light curve SNe la found pref- 
erentially in lower SSFR and/or more massive host g alaxies. Such effects must be accounted for 
when using SNe la as cosmological probes; for example. ISullivan et al.l (|2010l ) suggests the inclusion 
of an additional parameter in cosmological analyses to remove the host dependence. 



5.3. IMF Evolution 

By comparing the ages of the stellar population in the host galaxies with the age of the 
Universe at the redshift of the SNe, we can identify the first epoch at which star-formation occurred 
in the host galaxies. This provides an upper limit to the formation of stars which might be the 
progenitors of the Type la SNe. Since it is generally thought that < 8 M stars are the progenitors 
of Type la SNe we have used the host galaxy SED analysis to show that these low mass stars 
were in place 3 Gyrs after the Big Bang and possibly as early as z ~ 5 albeit with significant 
uncertainties tha t are related to the un certainties associated with measuring stellar population 



ages (Figure IT0|) . iTumlinson et al.l (|2004h argue that the nucleosynthesis yields as estimated from 
the metal abundances in halo stars and the electron scattering optical depths from the cosmic 
microwave background are well matched by requiring an IMF at z ~ 6 which is truncated at 
10 — 20 M & rather than the requirement that the first stars are very massive (> 14OM ). However, 
if our results are confirmed and low mass stars are found at redshifts as high as 5, this would rule 



-25- 



out such a truncated IMF. 



6. Conclusion 



We have studied the host galaxies of a sampl e of 22 Type la Super n ovae (SNela) at z > 0.9 5 



from lGilliland et al.1 <|l99flh : Blakeslee et alJ (|2003 ); lStrolger et al.1 (j2004l ): Ifiiess et al. 



(2004 



20071 ). 



We use the broadband photometry from HST ACS BViz, Spitzer IRAC as well as UJHK ground- 
based and some HST NICMOS JH data from the GOODS survey. We fit the photometry to the 
single stellar population models of Chari ot fc Bruzual (priv. com .) which are generated from 
the latest version of the GALAXEV code ([Bruzual Charlotl \2003i ) . We use Spitzer MIPS 24/rai 
data to place upper limits on the far-infrared luminosity of the hosts to break the well-known age- 
extinction degeneracy associated with SED fitting. We find the best-fit model for each host using 
a minimum x 2 technique to estimate the age of the stellar population of the SN la progenitors and 
hence place upper limits on the possible SN la delay times. We find evidence for both prompt and 
delayed SNela. When the distribution of stellar ages in field galaxies is factored in as a control 
sample, we find that the prompt SNe are dominant over the delayed SNe and that the SN la delay 
times possibly have a bi-modal distribution with a paucity of SNe with delay times of ~0.8 Gyr. 
As a consequence of the stellar population ages, we also show that the < 8M SNIa progenitor 
stars are in place by z ~ 2 and possibly by z ~ 5 (although with significant uncertainty) arguing 
against a truncated IMF in the first Gyr after the Big Bang. 
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Fig. 4. — Contour plots of the x 2 distribution for two representative SNela host galaxies (top: 
1997ff, bottom: 2003aj). The contours were calculated by collapsing over the remaining axes, 
i.e. for a fixed pair of age-Ay values the minimum x 2 allowing all other parameters to vary was 
found. This was then binned together. The green area is the 68% confidence level and the red 
area is the 95% confidence level. The confidence levels were found by finding those bins with 
X 2 < Xmin + X 2 gaussi where Xmin ls ^ ne minimum x 2 across the whole parameter range (i.e. the x 2 
of the best-fit) and Xgauss 1S the X 2 value which gives the 68% or 95% confidence level for a normal 
X 2 distribution with the same number of degrees of freedom as used in the fit. The blue cross gives 
the position of the best fit. 



- 32 - 



HST04P It 



0.95 < z < 1.2 
1.2 < z < 1.4 
1.4 < z < 

- 1.6 < z < 

2002k! 



1.6 

1.8 



X 



2 - 







1E-3 



HSTCWEog 



19971a 



S^ro' 



TO5Fer 
"HST05Lon 



HSTO40mb 
2005oj 2005dy 



200. iok ^STO4Mcg 



2002hp 



1E-2 0.1 1 

Luminosity Weighted Age [Gyrs] 



10 



Fig. 5. — The plot shows the luminosity- weighted ages of the stellar populations in the hosts with 
the associated xt- There are some SNela which originate from very young hosts (< O.lGyr). The 
vertical line is the weighted mean age of 0.07 Gyrs. 
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Fig. 6. — The plot shows figure [5] but when using the age upper-limits rather than luminosity- 
weighted ages. The vertical line is the weighted mean age of 0.08 Gyrs. 
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Fig. 7. — A histogram of the luminosity weighted ages of the stellar population in z > 0.95 Type 
la SN host galaxies. The ages appear to show a bi-modal distribution with the younger population 
having a median age of 0.11 Gyrs and the older population having a median age of 1.9 Gyrs. The 
black solid line shows the distribution for all SNe and the red dashed line shows the distribution 
when only high confidence SNe la are included. 
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Fig. 8. — Histogram of the luminosity weighted ages of the field galaxy sample in the same redshift 
range as the SNe hosts, used in the DTD calculation. 
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Fig. 9. — Type la Supernova Delay Time Distribution (DTD) for z > 1 SNe. The black points are 
those obtained for the whole sample, the red points are that obtained for the high confidence SNe 
only. Also shown is the DTD obtained bv lTotani et al.l (|2008l ) in green. Although, the uncertainties 
are dominated by the small number of candidates, if the DTD were a smooth power-law, we would 
have expected to find at least two additional SNe hosts with delay times in the 0.2-0.6 GYr range. 



- 36 - 



X 



5.5 
— i — 



3.1 
— i — 



0.95 < z < 1.2 
1.2 < z < 1.4 
1.4 < z < 1.6 

1.6 < z < 1.8 



2 - 



1.6 
— i — 



1.2 

— i — 



0.9 



HSTtMP* . 



HST04Eog 



2003es _ HST04Th<f s l 05K ° e 



nST05Fer 
HSTMLon 



HSTCBStr • 

!-• 

HS1 MOmb 



2003oj 
2003dy 



2002hp 
I 



~20A3ok HST04MC9 
I l_ 



9i 7fq- 



3 

T sf [Gyr] 



Fig. 10. — The plot shows the first epoch of low mass star formation for each host with the 
associated xl- The vertical line is the weighted mean epoch of low mass star formation of 5.6 Gyrs. 
The < 8M Q progenitor stars of SNela are certainly in place by z ~ 2 and possibly by z ~ 5. 
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